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In this note we will discuss the energy dependence of particle ratio fluctuations in heavy ion 
collisions. We study how the inherent multiplicity dependence of ratio fluctuations is reflected in 
, the excitation function of the dynamical fluctuations. Specifically, we will show that the observed 

■ excitation function of dynamical K/tt fluctuations is consistent with the expected dependence on 

' the number of accepted pions and kaons in both the STAR and NA49 experiments. 
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INTRODUCTION 



Over the last severalyears event-by-event fluctuations of many observables have been studied in relativistic heavy 
ion collisions [l|, l2J, 13, U, la, la, lD, la, 13, III], llll, ll3j ■ The measurement of these fluctuations may reveal separate 
event classes, fluctuations associated with phase-transitions etc. (for a recent review see [HI). The first measurement 
of event-by-event hadron ratio fluctuations has been carried out by the NA49 collaboration, which analyzed the 
fluctuations of the kaon to pion (K/n) ratio in lead- lead collisions at a center of mass energy of ^/s = 17.3 GeV [H]. 
This measurement found surprisingly little fluctuations, of the order of a few percent, which could be explained by a 
combination of Bose-Einstein correlations together with resonance decays [H, [Hi . 

Subsequently, the N A49 collaboration has measured event-by-event ratio fluctuations for several values of the center 
of mass energy, ranging from — 6.3 GeV to ^/s = 17.3 GeV ||. In addition, at the Relativistic Heavy Ion Collider 
(RHIC) the STAR, PHENIX, and PHOBOS collaborations have measured fluctuations up to center of mass energies 
of t/s — 200 GeV so that excitation functions of many fluctuation observables are now available over a wide range of 
energies. While most excitation functions show only little energy dependence, that for the kaon to pion ratio exhibits 
a steep increase towards the lower energies. So far this increase could not be reproduced in either thermal model 
calculations [l3| nor with the microscopic transport model UrQMD [3| and thus has sparked quite some interest and 
speculations concerning the QCD phase transition [l8j]. Another calculation using the HSD event-generator [llj can 
describe the general increase of the fluctuations towards lower energies but fails to reproduce the very steep rise 
exhibited in the NA49 data. 

In this note, we will point out that the observable generally used in the discussion of dynamical fluctuations, 

. t" = rr~ — t~ ( 1 1 

' dynamical data ^mixed events V / 

has an inherent dependence on the multiplicity of particles that are used in the experimental analysis - these are 
the particles located within the phase space domain covered by experimental acceptance with particle identification 
capability. This has some nontrivial consequences for excitation functions extracted with fixed target experiments, 
such as NA49, where the acceptance changes considerably with beam energy. Specifically, in case of K/i: - ratio 
fluctuations, the dynamical fluctuations depend on the inverse of the accepted number of pions and kaons, and since 
their number decreases with beam energy this may lead to non-negligible corrections, as we shall show. Indeed in [3| a 
significant centrality dependence of these ratio fluctuations was found at top RHIC energies, consistent with an inverse 
multiplicity scaling. In addition for a fixed target experiment the actual acceptance changes with beam energy. Since 
it is the multiplicity of the accepted particles which determines the fluctuations, this needs to be taken into account 
if one studies an excitation function. 
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It is the purpose of this note to discuss the multiplicity dependence of particle ratio fluctuations and in particular 
those of the kaon-to-pion ratio. After a brief review of the underlying formalism governing fluctuations of (particle) 
ratios, we will discuss several ways to remove the inherent multiplicity dependence in the definition of the dynamical 
fluctuations, c| ynamical , Eq. [TJ We then apply them to the fluctuations of the kaon-to-pion ratio and discuss the 
energy dependence of this observable in this context. 



II. FLUCTUATIONS OF PARTICLE RATIOS 



Following fla| the fluctuations of a particle ratio A/B are given to leading order by 



<4>' 



with the definition 



Since 



{A ) 2 + (B) 2 2 (A)(B)) +0{d) 



5 A = A - (A) 
SB = B - (B) 



(2) 



Eq. [2] can also be written as 



(SA 2 ) = (A 2 ) -(A) 2 

(SB 2 ) = (B 2 ) {B) 2 
{5 A SB) = (AB)- (A) (B) 



(A 2 ) , (B 2 ) 2 _(AB) 



\(A) 2 (B) 2 "{A){B) : 
In absence of any correlation, ($A 2 } = (A) and (5B 2 ^ = (£?), and the scaled variance (Eq. [2]) reduces to 
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The difference between the scaled variance a 2 and the uncorrelated scaled variance Cu ncorre i at ed ^ s commonly referred 



to as the dynamical fluctuations, "dynamical' 



2 2 _ 2 
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where ^dynamical is the variable usually used by the STAR collaboration. Obviously in the absence of any correlations, 
^dynamical = 0, by construction. Introducing the scaled correlations 



Cab 



(SASB) - S AB (A) 



vOT) 



(5) 



the dynamical fluctuations, "dynamical' can ^ e written as 



dynamical - + J^C BB 



Cab 



(6) 
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We note that the scaled correlations, Cab, typically do not or only weakly depend on the multiplicity as we will show 
explicitly in the context of a resonance gas. Consequently, Eq. [6] already suggests, that a simple scaling of cr 2 ynamical 
with number of charged particles may not be sufficient. De pending on which of the scaled correlations dominates, 
^dynamical ma y either scale with 1/ (A), 1/ (B) or with (A) (B) or some combination of those. Furthermore, if 

the particle abundances differ considerably, say jgy -C 1, as it is the case for the kaon to pion ratio, the dynamical 
fluctuations will be dominated by the least abundant particle, even if the scaled correlations are of the same magnitude. 
This follows directly from Eq. El In this case a scaling with 1/ (A) should work best. 
Finally, quantum statistics gives rise to additional correlations [H, H3, HH 

with 

where (+) stands for Bosons and (— ) for Fermions and ua{p) dp is the number of particles of type A in momentum 
bin (p,p + dp). The correction term due to quantum statistics is typically of the order of a 5 — 10% for systems 
of consideration, resulting in 0(1 — 2%) effects for the dynamical fluctuations [IR [20|. After these general remarks 
about fluctuations of particle ratios and their scaling with multiplicity let us turn to the specific case of kaon-to-pion 
ratio fluctuations. 

III. MULTIPLICITY SCALING OF K/tv FLUCTUATIONS 

Let us now turn to the specific case of K/tt fluctuations. In this case the scaled variance, Eq. [21 is given by 



(O^f) _ f( {S Kf) f ((Sn) 2 ) (SKSn) 



a K/7T ~ (f)2 {K)2 + (?r)2 2 {K){w) j (8) 

and the dynamical fluctuations, Eq. El are 

dynamical = ^ + J^f™ ~ ^)J^ C 

Since K = K + + and 7r — ir + + 7r_ 

(K) = (K + +K_) 

(tt) = (n + +7r_) (9) 
({5Kf) = ((8K+Y) + ((8K_f) + 2(8K + 8K.) (10) 
{(Sir) 2 ) = {(5n + y) + {(5n_f) + 2(Sn + 5n_) 

(SKStt) = (8K + 8ir+) + (8K_8tt + ) + (8K + 8ir_) + (8K_8ir_) (11) 
so that the scaled correlations, Eq. are given by 

((SK+) 2 ) + ((6K-) 2 ) + 2 (5K+6K-) - (K+ + KJ) 



Ckk 

Cmr 



((Stt+) 2 ) + ((Stt-) 2 ) + 2 (Stt+Stt-) - (tt + + tt_) 



7T4- + 7T_ 



(SK + Stt + ) + (8K_8n + ) + (8K + 8tt_) + (8K_8tt_) 

^Ktt = , (12) 

y/(K + +K-)(n + +K-) 

We note, that the "diagonal" scaled correlations, Ckk and C nn , also contain cross-correlations between the positively 
and negatively charged kaons and pions respectively. Therefore, correlations introduced by resonances such as the 
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0-meson and the po-nieson will enhance the "diagonal" scaled correlations and thus the dynamical fluctuations. 
Resonances decaying into a kaon and a pion, such as the i^Q-meson will contribute to the off-diagonal scaled correlation, 
Cktt and will reduce the dynamical fluctuations. This will be different if charge specific ratios such as are 
considered. In this case only resonances which decay in either two i^ + -mesons or two 7r_ -mesons will contribute 
to the diagonal terms, while the ifg-mesons do contribute to the off-diagonal scaled correlation. Consequently, the 
fluctuations of the charge specific ratio is not necessarily the same as that of the ratio of sums of negative and positive 
kaons/pions. This is also seen in the data by the STAR collaboration 

It may be instructive to discuss the above equations in the context of a simple model system which contains n 7r ± 
charged pions and n K ± charged kaons as well as n po neutral rho-mesons, n w omega-mesons, phi-mesons, and nx* 
neutral Kq and n R * J^o- mesons - Via their decay channels these resonances will give rise to the various correlation 
terms in Eq.[TU In reality, of course, there are many other resonances contributing to a and to the scaled correlations. 
The ones chosen here are the lightest and most abundant ones which lead to correlations and thus should provide 
a rough idea on how the different terms contribute [13] ■ A detailed investigation would involve a full study in the 
hadron resonance gas model, which is not the purpose of this paper. In addition to correlations due to resonances 
there are effects due to quantum statistics (20L l22|. as already discussed. Here, for simplicity we work with classical 
(Boltzmann) statistics. Noting that the branching ratio BR((f> — > K + + K_) ~ i and BR (Kq — > K + + 7r_ ) = 
BR [Kq — > K- + tt + ) ~ I the average particle numbers in our simple model are given by 

2 / 

(K) = (K+ + KJ) = n K+ + n K _ + + - \ n K * + rig. 



2 / 

= (71+ + 7r_) = n 7r+ + n„_ + 2n po + 2n u + - [n K * + n R , 



and we obtain the following expression [H| for the (co) -variances in Eq. [TT] 

1 2 

<(^+) 2 ) = n K+ + + 2 nK 5 
1 2 

({SK^) 2 ) = n K _ + -nj + -n R , 



<(^ + ) 2 ) 



pa 



3 ni 



((57r_) 2 ) = n lx _ + n po +11^ + -nx* 

(SK+6K.) = X -n^ 

(5tt + Stt^) = n po + 
2 

(5K + 5ir_) = -n K * 
2 



(5K_8tt + ) 



leading to (see Eq. [TT]) 
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The corresponding scaled correlations, Eq. [5] are 



((SK) 2 ) = (5Kl + 5K 2 _ + 2(5K + 5K-)) = (K)+n 4> 
({Sir) 2 ) = (8-K 2 + + 8-K 2 _+2{8-K + 8-K-)) = {n)+2{n PQ +n u ) 

(SKSn) = I (n K * +n R ,^j 



Ckk 



(K) 

2 (n po + n^j) 



2 riK * + ne-* 

c K7r = - ; K ° (13) 
3 

We find, that the scaled correlations in our simple model indeed depend only weakly (if at all) on the multiplicity, since 
both the number of resonances and the number of kaons and pions are expected to scale roughly with the multiplicity 
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or volume of the system. In a thermal system at fixed temperature they would be constant. Furthermore, as already 
discussed the correlations introduced via the resonances affect all three scaled correlations. While the i^-mesons 
control the off-diagonal correlation term, (SKSir), both the </>- meson and the po and u) contribute to the diagonal 
parts. The former reduce the dynamical fluctuations while the latter increase them. Putting everything together, the 
scaled variance, Eq. [2j is given by 

2 

K 1 7T 

leading to 





dynamical - | ^2 + ^2 (jjQ ^ 

or in terms of the scaled correlations 

dynamical = ( Tw\ CkK + ~U\ C ™ ~ ^ JjW) (tt) ^ J ^ 

Evaluating scaled correlations for our simple model, Eq.Qljl for a temperature of T — 170 MeV and vanishing chemical 
potential, we get 

Ckk = 0.1 

= 0.36 
C K , = 0.13 

Obviously, in our simple model all the scaled correlations are of the same order of magnitude. In addition, adding 
more resonances will likely reduce C„„ as there are many resonances decaying into only one pion, which add to the 
denominator, (ir), but not the numerator of C n „. While it would be worthwhile to study these scaled correlation 
coefficients in a full hadron gas model, here instead we want to concentrate on simple phenomenological scaling rules, 
with a special emphasis on the effect of a varying acceptance. 



IV. PHENOMENOLOGICAL SCALING 



In this section we want to discuss several ways to scale out the multiplicity dependence of the dynamical fluctuations. 
Of course, if all the relevant scaled correlations are known and if they, as we argued, depend only weakly on the 
multiplicity and beam energy, the appropriate scaling is simply given by Eq. [TH This is equivalent to having full 
understanding of all the sources for the fluctuations in which case this discussion is mute. In general, however, we do 
not have a full understanding of all the sources contributing to the fluctuations. In this case the multiplicity /energy 
dependence may provide additional information such as signals for a possible phase transition etc. In order to extract 
this information, one needs to understand "trivial" dependencies on the multiplicity, as exhibited e.g. in Eq. [HI 
In the following we will discuss several "trivial" scaling prescriptions which we believe should be applied before 
conclusions about new phenomena can be drawn. We will focus on the rather interesting energy dependence of the 
JsT/7r-fluctuations, referring to 0] for a discussion on the centrality dependence. In addition to applying appropriate 
scaling prescriptions to the dynamical fluctuations, it is essential to realize that the multiplicities which control the 
fluctuations, such as (K) and (tt) are that of the identified particles and not the extrapolated total multiplicities or 
mid-rapidity multiplicities. This is especially important when studying the energy dependence of fixed target data, 
as the acceptance varies with beam energy. We will discuss this point in detail below after we have introduced the 
various scaling prescriptions. 



A. Multiplicity Scaling Prescriptions 

As suggested in [la ], one way to avoid any scaling of the dynamical fluctuations with the number of accepted 
particles would be to take the ratio of the measured scaled variance a 2 over that of mixed events cr mixcd , 

/--£-. 

(7 

mixed 
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instead of the difference as it is done in the definition of ^dynamical! Eq. [TJ This, however, has the disadvantage, that 
correlations and fluctuations due to the the detector do not cancel out. In order to remove the multiplicity dependence 
in the same fashion, one can simply divide the dynamical fluctuations, Od ynamical , by that of uncorrelated particles, 
^uncorrected; Eq. [31 evaluated for the same number of particles. 

2 

t - ^dynamical (15) 



J oisson 



a 2 



uncorrelated 



with 



'"uncorrelated / t*-\ 



(K) (tt) 

for the case at hand. This scaling we will subsequently denote as "Poisson" scaling, as we scale with the scaled variance 
of a Poisson distribution based on the observed multiplicities. The advantage of this scaling is that it is unbiased in 
the sense that one does not need to make any assumptions about the relative magnitude of the scaled correlations. 
Another unbiased scaling would be to scale with the number of particles involved, i.e. N p = (K) + (tt) . This is similar 
in spirit of 0] where a scaling with the number of charged particles was studied. We shall henceforth refer to this as 
particle number scaling, 

/Particle Number = {(K) + (tt)) CTdynamical (16) 

Finally the expression for cr^ ynamical in terms of the scaled correlations, Eq. [HI suggests the scaling with the kaon or 
pion number or with the geometric mean of both, depending on which of the scaled correlations dominates. 



/k aon Number — 

(K) 

dynamical 
dynamical 



/Pion Number — (^) <7 d 



/geometric = \J (K) (n) 

^"dynamical 

Since the number of kaons is much smaller than the number of pions, at least for the lower energies the kaon- number 
scaling is similar to the Poisson scaling, 

1 1 _ (K) + (tt) _ 1 

^correlated ^ + ^ ^ ^ - ^ . 

for (K) <C (tt) . Alternatively the above scaling relations allow to relate the dynamical fluctuations at a given center 
of mass energy with those at another energy. Specifically, given the dynamical fluctuations at top RHIC energies, 
^fs = 200 GeV, the dynamical fluctuations at any other center of mass energy is given by 



Poisson scaling: 



I r\ ^nnr, ,rt \/ (K) + (I) 
^dynamical [V s ) = f dynamical (2UU(jeVJ — = (17) 



(K) ' (tt) 

• Particle Number scaling 



i i l 

200 GeV 



t tt rt nnP„,A #) + \7T/|200GeV „ s 
^dynamical (VSJ = CTdynamical (^00 UeV j =r (18) 



Nk -scaling: 

^dynamical (Vs) = f dynamical (200 GeV) V N j^=^p ^- (19) 




N„ -scaling: 



I r\ toaari \T\ V^Tl 200 GeV , . 

f dynamical (VSj = CTdynamical (ZULUjeVj -= (20) 
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• Geometric Scaling: 

f dynamical (V*) = ^dynamical (200 GeV) — - — (21) 

The resulting scaled dynamical fluctuations are shown in Fig. Q] where we find that the energy dependence of o"dy™mical 
is reasonably reproduced by any of the scaling rules discussed above. The values for the multiplicities of identified 
particles taken from 0,0] as well as the result for the rescaled fluctuations according to Eqs. i|17H21|l . The essential 
point for the success of these scaling rules is that we have used the number of identified kaons and pions for the 
mean values, (K) and (ir), entering the scaling formulae, Eqs. [T71I2T1 This leads to an additional energy dependence 
especially for the NA49 data. Since NA49 is a fixed target experiment, the actual acceptance and thus the fraction of 
identified particles of the total number of particles may vary considerably with the beam energy. This is illustrated 
in Fig. [2j where we applied the same scaling formulae, Eqs. [171I21) but used the mid-rapidity yields, = 0) for 

the respective mean particle numbers, (K) and (71") [13, 0, Obviously the energy dependence is not reproduced, 
especially for the highest SPS energies, -y/s = 12.3 and 17.3 GeV. 
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100 



E (GeV) 

cm v / 



Figure 1: Different scaling scenarios, Eqs. [1711211 based on the y^s = 200 GeV data from STAR 



Looking more closely at Fig.[H we see that the STAR data at y/s — 62.4 GeV, which have a rather small error-bar 
are not well reproduced. Instead of arguing for new physics in this energy regime, we note that in Ref. the two 
most central values for <Xd yn amicai at this energy do not agree very well with the systematics developed by the STAR 
collaboration either. As already discussed in section HU since (K) / (tt) <C 1 the dynamical fluctuations are dominated 
by the kaons. Therefore, the kaon-number or Poisson scaling should work better than pion-number or particle number 
scaling. This is indeed the case. 

One may be tempted to use the quality of agreement of the various scaling prescriptions to draw conclusions about 
the importance of the various contributions to the dynamical fluctuations. Given the experimental error-bars and 
the quality of agreement of the various scaling rules the value of such an exercise is not obvious. At least, a global 
fit based on the various scaling prescriptions needs to be carried out before any more detailed conclusions about the 
strength of the various contributions can be drawn. 

Finally let us comment on the fact that calculations based on the UrQMD model do not reproduce the rise of 
the dynamical fluctuations as observed by NA49, although the NA49 acceptance has been applied First of all, 
while UrQMD seems to do reasonably well for mid-rapidity abundances, it is not clear if UrQMD does reproduce the 
correct multiplicity of identified particles within the acceptance, especially at small center of mass energies. Second, 
it is not at all obvious, if UrQMD contains indeed all the relevant sources for fluctuations and correlations. For 



8 





E (GeV) 

cm v ' 



Figure 2: Scaling with -^{y = 0) instead of the actual number of identified particles contributing to (Xdynamical- 
same scaling formulae are used as in Fig. [T] 



Otherwise 



instance quantum statistics is not taken into account. The fact that the above scaling relations are able to connect 
the dynamical fluctuations over a wide range of beam energies does not say anything about the nature and origin of 
these fluctuations except that they seem to be the same at all energies. While this considerably weakens the argument 
for the isT/7r-fluctuations being a signature for the QCD critical point, the disagreement with UrQMD may very well 
point to a yet to be discovered new source for fluctuations/correlations. This however, requires that all "trivial" 
sources, such as e.g. quantum statistics, are systematically taken into account. 

Finally, let us conclude this section by proposing an observable, which should, to leading order, be independent of 
the multiplicity. As already discussed at the beginning of this section, the ratio 



mixed 

proposed in [l5| has the disadvantage that correlations and fluctuations due to the detector do not cancel out. We, 
therefore, propose to study instead the ratio, Eq. [TBI 

2 

p ^dynamical 

JPoisson — 2 

^uncorrelated 

or equivalently /p isson + 1, which has the same limit as / in the absence of correlations. In Fig. [3] we show /p isson + 1 
for both the K/ir fluctuations as well as for the p/7r-fluctuations, which have been measured by NA49 as well. We 
see a rather weak energy dependence of the if/-7r-fluctuations, whereas the p/7r-fluctuations exhibit a variation with 
energy. We note, however, that the energy dependence of the p/7r-fluctuations in the scaled variable is opposite to 
that in "dynamical- While 0d ynam j ca i decreases with energy, the magnitude of /p isson increases, suggesting in increase 
of the the strength of the correlations with energy, which seems to level off at top SPS energies. This suggests that 
the correlations, mostly due to baryon resonances, increase with energy up to y/s ~ 15 GeV. Using the chemical 
freeze-out parameters of Ref. [IB] and assuming that only the delta resonance contributes to the scaled correlation 
coefficient C p7r we find indeed an increase of the correlations from AGS up to SPS energies. 



V. CONCLUSIONS 



In this article we have reviewed the multiplicity dependence of particle ratio fluctuations. We have provided 
several scaling prescriptions which correct for the inherent dependence of the dynamical fluctuations "dynamical on 
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Figure 3: The rescaled fluctuations /poisson + 1 for K/tt (black) and p/-7r-fluctuations (red). 
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i 71 ) ident. 
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0"dyn.[%] 


Geom. scaling 


TV^-scaling 


-scaling 


Poisson-scaling 


Part. Num. scaling 


/poisson(%) 


6.27 


30.9 


5.7 


7.89 


10.05 


7.15 


14.11 


7.61 


13.27 


3.0 


7.63 


66.6 


10.2 


6.90 


7.17 


5.35 


9.60 


5.61 


9.15 


4.2 


8.77 


103.3 


15.0 


5.47 


5.83 


4.41 


7.71 


4.61 


7.37 


3.9 


12.3 


227.6 


31.2 


3.70 


3.99 


3.05 


5.19 


3.18 


4.98 


3.8 


17.3 


416.6 


54.2 


3.17 


2.99 


2.32 


3.84 


2.41 


3.69 


4.8 


19.6 


227.6 


10.7 


4.04 


5.21 


5.21 


5.19 


5.21 


5.19 


1.7 


62.4 


319.2 


14.5 


3.61 


4.44 


4.48 


4.38 


4.48 


4.39 


1.8 


130 


351.9 


14.4 


3.78 


4.34 


4.51 


4.18 


4.49 


4.19 


2.0 


200 


432.6 


20.6 


3.77 


3.77 


3.77 


3.77 


3.77 


3.77 


2.8 



Table I: Table of the results for the various scaling scenarios Eqs. I|17II21I) using the STAR data at ^fs 
The data for identified multiplicities, (7r) ident and {K) ident , are from [3|, [7|. 



200 GeV as reference. 



the number of identified particles. We have demonstrated that these scaling rules naturally reproduce the trend seen 
in the energy dependence of the kaon-to-pion fluctuations. Consequently any interpretation of the rise of adynamics! 
towards lower energies in the context of a possible QCD critical point, needs to account for the "trivial" effect due 
to the multiplicity dependence. We propose that future studies of energy and / or multiplicity dependencies should 
correct for the "trivial" multiplicity dependencies inherent in Cdynamicai- In our view, the simplest and least biased 
scaling is the Poisson scaling, /p i SS on, which could serve as a benchmark. Of course a more complete analysis of the 
energy dependence would include a thorough global fit of the available data, which we have not carried out in this 
paper. 

We have further applied the scaling to the p/7r-fluctuations measured by the NA49 collaboration. We find that the 
properly scaled observable still exhibits a strong energy dependence, which, however, is opposite to the unsealed one. 
In this context it would be interesting to include the STAR data, which unfortunately are not available in a final 
version yet. 

Finally we point out that any multiplicity scaling needs to be based on the mean multiplicity of the actual identified 
particles used for the fluctuation measurement instead of an extrapolated multiplicity such as AN /Ay. 
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